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fH , We provide a complete and consistent framework to include subasymptotic 

Q-i! perturbative as well as mass corrections to the leading twist (r = 2) evalua- 

Q^l tion of charged and neutral current weak structure functions and the pertur- 

bative neutrino cross sections. We revisit previous calculations in a modern 
language and fill in the gaps that we find missing for a complete and ready- 
to-use "NLO ^-scaling" formulary. In particular, as a new result we formulate 
the mixing of the partonic and hadronic structure function tensor basis in the 
I operator approach to deep inelastic scattering. As an underlying framework 

we follow the operator product expansion a la Georgi & Politzer that allows 
the inclusion of target mass corrections at arbitrary order in QCD and we 
provide explicit analytical and numerical results at NLO. We compare this 
approach with a simpler collinear parton model approach to ^-scaling. Along 
with target mass corrections we include heavy quark mass effects as a cal- 
culable leading twist power suppressed correction. The complete corrections 
have been implemented into a Monte Carlo integration program to evaluate 
structure functions and/or integrated cross sections. As applications, we com- 
pare the operator approach with the collinear approximation numerically and 
we investigate the NLO and mass corrections to observables that are related 
to the extraction of the weak mixing angle from a Paschos-Wolfenstein-like 
relation in neutrino-iron scattering. We expect that the interpretation of neu- 
trino scattering events in terms of oscillation physics and electroweak precision 
physics will benefit from our results. 



I. INTRODUCTION 



Recent neutrino experiments have shown strong evidence for neutrino masses and mixing. 
A variety of experiments [1,2] have yielded data that when combined [3,4] have led to a 
picture of neutrino mass and mixing with large mixings between <-> i/^ and v^j, 
Am^2 ~ 7 X 10~^ eV^ and Am^g ~ 2 x 10"^ eV^. Solar neutrinos are produced with E^, < 20 
MeV [5], making z/g ly^ the relevant oscillation process. For the atmospheric case, (i?^) ~ 1 
GeV and — > i^r is the dominant oscillation process. The MINOS, OPERA and Chorus 
experiments will study z/^ Uj. oscillations, either by disappearance or Ut- appearance 
[6]. 

In the GeV neutrino energy range, target mass effects are important. For muon neutrino 
charged current interactions, at low energies {E,, ~ 100 — 800 MeV), the quasi-elastic process 
dominates and nucleon mass effects are easy to incorporate [7]. At intermediate energies, 
exclusive few pion production processes are the largest contributions to the cross section 
[8-10]. Above a few GeV, deep inelastic scattering (DIS) dominates the cross section. The 
energies for tau neutrino charged current interactions are higher because of the tau lepton 
threshold. A scries of conferences (NUINT) has grown out of the efforts to improve and 
combine knowledge about the distinct neutrino interaction modes into event generation for 
oscillation searches [11]. 

In a different branch of recent research at the interface of neutrino physics and QCD, the 
NuTeV collaboration has measured the weak mixing angle from charged and neutral current 
muon neutrino and anti-neutrino interactions with a broad band beam (see Fig. 5). The 
beam flux maximum is at an energy of about E^ ~ 60 GeV which translates into an average 
event energy of about E^, ~ 100 GeV because the cross section rises approximately hnearly 
with energy. NuTeV extracts a value of sin^ Ow = 0.2277 ± 0.0013 ± 0.0009 [12] which lies 
about 3(7 above the standard model value 0.2227 ± 0.00037 [13]. The latter fit value does 
not include neutrino scattering data. Possible explanations of the NuTeV anomaly within 
and beyond the standard model have been mapped out in Ref . [14] . Nuclear physics related 
complications in its interpretation - the NuTeV target material consists mostly of steel - 
have also been discussed in Ref. [15]. No consistent picture has as yet emerged from these 
efforts to understand the measured value. 

The perturbative QCD parton picture is applicable at NuTeV energies and at the higher 
energy end of the oscillation neutrino beams. A perturbative background to (quasi-) elastic 
scattering is present at lower energies as well. At none of these energies can sub-asymptotic 
corrections to the simplest parton model picture be neglected, and both oscillation and 
precision neutrino physics arc bound to benefit from a concise inclusion of these terms into 
data analysis. At least for precision physics, it does not suffice to model these corrections but 
it is a necessity to extract theory parameters using theoretically sound cross sections. E.g. a 
well defined separation of higher twist effects (sec e.g. [16]) requires that power suppressed 
mass terms are accounted for correctly in the leading twist-2 cross section. 

In the present article we include 

• 0{as) perturbative QCD corrections 

• O(M"/Q'0 target mass effects 

• 0{mP' IQ'^) and 0{\n.m? /Q"^) heavy quark mass effects 
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• 0{'m'i"/{M'^E")) heavy lepton mass effects (mostly for r production) 

and any combination of the above. One could summarize the above as (Nachtmann or Georgi 
& Politzer) {-scaling for weak structure functions merged with NLO QCD for light and heavy 
quarks: While the heavy quark and perturbative corrections come in through the expansion 
of the Wilson coefficients in the operator product expansion, the target mass corrections 
find their way into the final result through the Lorentz structure of the corresponding (non- 
reduced) operators. 

Target mass corrections (TMC) to DIS have been know for a long time. The first 
discussion in terms of the operator product expansion (OPE) at leading order in QCD was 
done by Georgi and Politzer [17] in 1976. Later, these same target mass corrections were 
derived from a parton model point of view by Ellis, Furmanski and Petronzio [18]. DeRujula, 
Georgi and Pohtzer discussed next-to-leading order (NLO) QGD corrections to target mass 
corrected structure functions Wi^2 in the context of local duality in electroproduction [19] 
and using off-shell regularization. In this paper, we present the target mass corrections 
for charged current (CC) and weak neutral current (NC) and Ut DIS with nucleons, 
including next-to-leading order QCD corrections and heavy quark production using modern 
conventions. 

In the next Section, we outline the procedure for evaluating the TMC beyond leading 
order in QCD within the OPE. In Section III, we show our results in analytic form and 
compare them to the collinear parton model. We also make a few comments on nuclear 
targets and on the validity of the Albright & Jarlskog relations for the structure functions 
F4_5. Numerical results for the structure functions are shown in Section IV for both the 
OPE and parton model approach to target mass corrections. As a further application, we 
evaluate observables that are related to the NuTeV Weinberg angle analysis in NLO QCD 
with high numerical precision. A summary appears in Section V. 



II. TARGET MASS CORRECTIONS 

A. Parton Model Approach (Briefly Revisited) 

Target mass corrections appear in a number of ways that can be described in parton 
model language. First, the parton fraction { (Nachtmann variable [20]) of the light cone 
momentum of the nucleon = (P^ + P^)/\/2 is related to the Bjorken x by 

1 1 n M2" ^ 2x , , 

^ (2-1) 



e 2x V4x2 g2 s 



for massless quarks. Second, there is a mixing between partonic structure functions in 
the evaluation of the hadronic structure functions. The partonic light cone momentum is 
related to the hadronic light cone momentum via p"*" = {Pj^, however, partonic momentum 
projections are not simple rescalings of the hadronic momentum projections on the hadronic 
tensor W^^ because p~ ^ ^P^ when {p^,P^} = {0,M^}. Finally, in a collinear expansion 
[18] target mass effects appear from transverse momentum (/c^) effects. 

In the collinear parton model (/c^ = approximation), target mass corrections have 
been evaluated in Ref. [21] and later in Ref. [22]. The intrinsic kx of the target partons 
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is limited by the nucleoli mass M, which introduces further corrections. Ellis, Furmanski 
and Petronzio [18] showed the equivalence of a non-coUinear parton approach results, where 
the parton is on-shell but not coUinear with the nucleon, to the operator product expansion 
target mass corrections, to which we will turn next. 



B. Operator Product Expansion 

The three contributions of the partonic approach appear automatically in the opera- 
tor product expansion approach to target mass corrections, first discussed in this context 
by Georgi and Politzer [17]. In this approach, one starts with the relation between the 
hadronic tensor W^" and the virtual forward Compton scattering amplitude T^'^. For weak 
interactions, the relevant set of structure functions is Wj, i — 1 — 5: 

W,.= ^ / e"^-^dh{N\ [J,(z), J,(0)] \N) (2.2) 
= - Disc / e'^'-'d^ziNliT {J^{z)J^{0)) \N) (2.3) 
= - Disc T^, (2.4) 

TT 

- -g^.W^ + -j^W2 - ie^,p,^W^ + -j^W^ + — . (2.5) 

In Eqs. (2.3), (2.4) the discontinuity is taken across the v = > Q positive frequency cut: 
Disc T^, = lim 1 [T^,{u + is) - T^,{v - ie)] . (2.6) 
The hadronic tensor W^y is probed in deep inelastic interactions with a leptonic current 

4 = ^2liJ.{'^i - a«75)V'i (2.7) 

through exchange of a boson of virtuality Q^. The corresponding cross section in the target 
rest frame is 

^ = Al + QVMlY EM^^y^K,v„a,rm,m,)FM ) (2-8) 

where Mb=w^,z° is the boson mass, x and y are the standard DIS variables and where the 
modern normalizations are 

{Fu F.. F„ F,, F.} = [Wu ^W,, ^W,, ^W., ^W,] . (2.9) 
For purely electromagnetic interactions* we have 



*Interference contributions between 7* and at very high (e.g. at HERA) are not sensitive 
to target mass corrections. They are, therefore, not covered by Eqs. (2.8), (2.10). 
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dx dy Q'^ 



Y,Wi{x,y,Eul,Q,mumi)F,{x,Q'') . (2.10) 



j=i 



The electroweak weights Wi in (2.8), (2.10) are hsted in Eqs. (Al)-(A15) for general lepton 
masses and couphngs along with the experimentally relevant cases. 

As is well known, contributes only to weak processes because of the parity violating 
nature of the anti-symmetric e-tensor. A sixth independent tensor combination ~ WQip^q^,— 
PuQi^) does not couple to the the most general leptonic current 

L'"'{vi,ai,mi,m2) iPuQu - Pi^Qn) = y{vi,ai,mi,m2} (2.11) 

where L^'^ is the leptonic analogue of the hadronic W^^/. Because We decouples from the 
DIS process, we omit it from the beginning and throughout. 

In neutral current reactions, W^^^ contribute through the axial coupling ai only, i.e. they 
do not couple to a 7* exchange. In standard formulae, the tensors corresponding to W4 
and W5 are sometimes combined with those corresponding to to Wi and W2 so that cur- 
rent conservation is manifest in the electromagnetic case. Generally, for both neutral and 
charged current weak reactions the structure functions W4 and W5 do not appear in the 
usual discussion since they are multiplied by the ratio of the lepton mass squared to ME^. 
They obey the Albright- Jarlskog relations [23] which are valid at leading order in the limit 
of massless quarks and a massless target^: 

F4= (2.12) 
2xF^= F2 (2.13) 

The forward scattering amplitude, explicitly including perturbative corrections, can be 
written as 

00 

T'''= E [-9'''%A,.Cf + g^.X.Q'C? - i^'"'''^9a,M,.Cf (2.14) 

k=l 



-g^gw^M^Cf + «9'g/.2 +<lg''^//.2)C5''■)g^3•••9mgife- 



where A2k represents the reduced (scalar) matrix element of a twist-2 (traceless and sym- 
metric) operator Qf^^::^'^^'' of spin 2k 

{N\ O'^li"^^ \N) = ^2itn'^'-'^2fc _ (2.15) 

The QCD perturbative corrections come in through non-trivial Wilson coefficients In 
Eq. (2.14), the factorization into universal operator matrix elements and process-dependent 
coefficient functions is a consequence of the operator product expansion-!- [24] of the product 



^The conditions under which Eqs. (2.12) and (2.13) are valid individually are detailed in the 
appendix of Ref . [22] and in Section III D 2 below. 

■i-The OPE has found its way into many excellent textbooks and clear pedagogical introductions 
can be found in e.g. [25]. 
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of currents T^^ in Eq. (2.3). We assume that O^l^^''^* of twist 2 to make contact with 
the parton model. For the following evaluation of the target mass corrections, however, it is 
sufficient that C is a traceless and symmetric operator. Symmetry under permutation and 
tracelessness under contraction of any pair of indices of O^Jf^^^^fc require then that the tensor 
part of the matrix element in Eq. (2.15) is 

n'^^---'^- = E(-i)^'^§^^__^ (/)^ (2.16) 

where g ... g p ... p abbreviates a sum over {2k)\/[2^ j\{2k — 2j)!] permutations (not count- 
ing identical terms such as g^'^^^ — g^^^^ twice) of the indices. In Eq. (2.16) the j — 
approximation Ht^^-'t^'^k oc p^i....p^2fc reproduces the massless parton model and the j > 
terms resum the target mass corrections. While they are formally power suppressed of order 
O^M"^^ / Q^^), the corrections derive from the same operators Qf^^::^'^'^'' as the j = term 
and as such are of twist-2. This is in contrast to power corrections from operators of higher 
twist C'^^^'^^* which introduce new non-perturbative matrix elements independent of those 
in Eq. (2.15). In classical terminology twist classifies strictly the dimension minus the spin 
of an operator; the difference between target mass corrections and higher twists is, neverthe- 
less, sometimes referred to as a distinction between kinematic and dynamical higher twists 
as well. 

The quantity is the 2k^^ moment of a universal function / or the (2 A; — 1)**^ moment 
of the quark distributions q{x)^ 

= C dx x^^fix, Q^) = C dx x^'-\{x, g2) . (2.17) 
Jo Jo 

Beyond leading order we include the non-universal Wilson coefficients into the moment 
integral 

C^'A2k= f'dxx^^ Mx,Q^) , (2.18) 
Jo 

thus defining a set of non-universal x-space functions /i=i,...5. The fi can be related to the 
zero mass limit of the experimental structure functions from the j — term of the 
series (2.16) as will be explained in the next section; the effect of the resummation of the 
J > terms will induce convolution integrals over the f/^-* and lead to Nachtmann scaling 
[20]. 

The standard calculation a la Georgi-Politzer of extracting the moments of the struc- 
ture functions, then performing inverse Mellin transforms, leads to equations most simply 
expressed in terms of /i(^, Q^) and 



^For ease of notation, throughout this article we sometimes suppress scale dependence as a 
functional argument. Q^-dependence enters through the renormalization of (Jf^]^;^^'"' i.e. strictly 
^2fe = ^2fe(Q^)- Also, for the purpose of this article and in order to make contact with the parton 
model it suffices to assume that O^^^'^^fc jg ^ quark operator since the inclusion of operators built 
from gluons and operator mixing would not alter the discussion anywhere. 
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(2.20) 



(2.19) 



The middle part of Eq. (2.20) is a more transparent representation when we write down the 
results below. For numerical evaluation, however, it is more convenient to express the double 
convolution integral as a moment of a single convolution as suggested by the right-hand-side 
of the same equation. 

We have extended the Gcorgi-Politzcr [17] and DeRujula-Georgi-Politzer [19] analyses 
by including the full set of weak structure functions along with their Wilson coefficients Cf^. 
Abelian gauge invariance is not assumed through a priori relations among C^^, although it is 
restored in the electromagnetic scattering case by the explicit values of C"^^. In unpolarized 
DIS, the leading component is of twist-2 for any Wi=i,...,5, and the coefficient functions do 
the bookkeeping of the tensor components in W^v One finds a mixing of structure functions, 
which for Fi was already found in Ref. [19]. The general pattern of mixing can be easily 
inferred from the j > terms in Eq. (2.16) that add different tensor components to the j = 
term. We note that in polarized DIS, the spin dependent structure functions [26-29] gi and 
g2 receive contributions from twist-2 and twist-2 and 3 operators, respectively. In this case, 
tracking the Wilson coefficient basis follows automatically from tracking the operators of 
twist two and three and from the Wandzura-Wilczek relation for the twist-2 part [30]. 

The method of evaluating target mass corrections as inverse Mellin transforms of mo- 
ments of structure functions was presented in explicit detail in the Georgi & Politzer analysis 
[17] already. Recent analogous calculations in the spin dependent case [28, 29] have repeated 
the mathematical steps involved. Apart from the differences listed above, we feel it is un- 
necessary to add a further derivation to the literature. In the next section, we show our 
explicit results for the target mass corrected structure functions Wi and the corresponding 



The procedure outlined above leads directly to the results for the structure functions in 
terms of Bjorken x, the Nachtmann variable ^, and as well as = hP jQ^ . The structure 
functions are: 



III. RESULTS 



A. Generic Formulae 




,2 Q \ 92M) 

(9x L 1 + 



(3.1) 



(3.2) 



(3.3) 



dx^ 1-1- //^^ 




ax2L(i-;,e)(i+/^e) 



(3.4) 
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W^— —fix" 



d_ 
dx 



2^3 



dx^^i-fiexi+f^e) 



After performing the derivatives in Eqs. (3.1)-(3.5), one gets a generic form 



(3.5) 



(3.6) 



=1,5 



The generic formula for the Wj apphcablc to Eq. (2.5) is more useful translated to 
modern normalization conventions of the structure functions (Fj) or to the parton model 
language. In the next sections, we associate the fi and related hi and g2 to the experimental 
structure functions Fi{x,Q^) in the M — > limit. The Fi{x,Q^)\j^^Q can be evaluated at 
NLO using parton distribution functions as documented in detail in [22,31]. As in [22] and 
as is appropriate for the neutrino energies under consideration, wc will be working in the 
MS factorization scheme with Uf = 3 active flavours. To avoid redundancy, we will not 
reproduce the formulae of [22,31] in the present article. In Section lllC we comment on 
nuclear target corrections. We then consider the association of Wj to parton distribution 
functions in Section HID. The master formulae for Fj in structure function language and 
in parton model language appear in Eqs. (3.17) and (3.27) below. 



B. Target Mass Corrected Structure Functions 

Our procedure for determining the quantities fi and associated hi and Qi is to consider 
the experimental structure functions Fi in Eq. (2.9) with the Wi in Eqs. (3.1)-(3.5) in the 
— > limit 

Fj'-^-* = , (3.7) 

in which ^ ^ Note that the Fj*-"-* here denote functional forms and not observables - 
though they are accessible in principle through the inversion of Eq. (3.17) below as the 
structure functions with moments of definite spin in the Nachtmann sense [20] or the per- 
turbativcly corrected scaling functions in the duality sense [19]. To make the distinction 
from more obvious we will, below, sometimes label the measurable structure functions 
Fj™*". From the j = term of the target mass series in Eq. (2.16) we find the following 
relations between fi and the zero-targct-mass limit functions F/°^ 

{Fi'\x),Fi'\x),Ft\x),Fi'\x),Ft\x)} = [pi{x),x'f2{x),xfs{x), ^f,{x), |/5(x)} . 

(3.8) 

In Section III.C, we show that this form explicitly substantiates the partonic interpretation 

fiix) = ki fix) + O {as) ^ki^ + O (as) (3.9) 
with the couplings of the quark current 

= ^2lix{vq - aq-tb)A (3.10) 

entering as 
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h = vl + al (i=l,2,5) 
ki = 2 Vqttg (i=3) 
ki = (i=4) . 

The form of Eq. (3.9) is modified slightly to include quark masses. For the applications in 
Section IV we calculate the Fj^'' using the 0{as) coefficient functions in Refs. [22,31]. Here 

we focus on relating F™'-'' to the fJ°^ and their convolutions. With Eq. (3.8), Eqs. (2.19) 
and (2.20) now become 

hi^,Q')=l'de^-^^^^^^^ (3.11) 
h2{^,Q')=l'de^^^^^^ (3.12) 
h{^,Q')= I' de ^^^^^ (3.13) 

h,{LQ')= di' ^-^^^^^ (3.14) 
h.{i,Q')= /^e'^^f^l^ (3.15) 



92{i, Q')= f de h,{i', Q') = f di'ii' - (3.16) 



Prom Eqs. (3.1)-(3.5), (3.8), (2.9), (3.11)-(3.16), the target mass corrected structure func- 
tions then are 

if MC ^ ^ ^jFf (e, Q') + B]h{^, Q') + CM^, Q") ■ (3-17) 

1=1,5 

The coefficients {A'j, Bj, Cj} are listed in Table I-III in terms of x, ^, and p where 



Note that 



1 



X = . (3.19) 



While we have sometimes referred to twist-2 operators to meet the parton model, Eq. (3.17) 
holds for contributions from any traceless and symmetric operators and it is independent of 
the perturbative order in Ug. For leading order 0{a^) twist-2 coefficient functions, Fi^2 as in 
Eq. (3.17) agree with the standard results in Ref. [17] and the non-coUinear parton results 
in Ref. [18]. For F^ there is a factor of 2x mismatch in relative normalization of the F^^^ 
and hs term between Ref. [17] and Ref. [18]; we agree with the results in Ref. [18]. 
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C. A Note on Nuclear Targets 



By the very nature of weak interactions heavy nuclei are more natural targets than a 
nucleon, so a brief comment on large Ma = AMjq nuclei may be worthwhile. In the OPE the 
coefficient functions are manifestly target independent and the target dependence resides in 
a generalization of Eq. (2.15) with lA^") replaced by a general nuclear target 

(Z, A\ C)^^-'^^. \z^A) = 4f''^)n^i-'*2fc _ ^3 20) 

i.e. ultimately the nuclear dependence for r = 2 is contained in nuclear modifications of the 
PDF in Eq. (2.17)** 

= A f'dx x^^-\^^^^\x, Q2) . (3.21) 

Genuine higher twist r > 2 effects may be of different relative size in nuclei, e.g. the twists 
investigated in Ref. [32] are enhanced by a factor A^^^. The M\ = A'^Mff enhancement 
of target mass terms is, however, balanced everywhere by kinematic factors that scale with 
inverse powers of A such as 

_ _ _x 

" 2P^ ~ 2APr, ■ q ~ A ^^'^^^ 

Neglecting for the moment the nuclear modifications in Eq. (3.21), under the assumption 
that^t 



^^''^\x) = A q{x) e{l - x) (3.23) 



target mass corrections are then the same for nuclear targets. Eq. (3.23) neglects the large 
1 < X < A tail in nuclei which is (small but) nonzero due to Fermi motion. 



D. Comparison with the Collinear Parton Model 

So far, no reference to the parton model was required in deriving the results. It is 
interesting to compare the above results, derived from the operator product expansion, with 
the collinear {kj_ = 0) parton model approach to target mass corrections [21,22]. The 
latter also leads to Nachtmann scaling but the mixing of the tensor basis comes in only 
from the non-collinear k" ^ component of the parton momentum. This approximation 
avoids single and double-convolution terms of the form of Eqs. (3.11)-(3.16). If the collinear 
approximation turns out to be of reasonable accuracy, one may hope to employ similarly 
approximate mass corrections to processes where the rigorous operator approach is not 
available such as processes where the final state is not fully inclusive. 



**Iii the absence of nuclear effects we could drop the label "Z, ^" on q(^'^^ in Eq. (3.21). 
ttThe isospin orientation is attributed to the PDF without explicit labeling. 
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1. The CC Process including Charm Production 



We now find relations between fi and the parton model structure functions for scat- 
tering off of the Cabibbo-Kobayashi-Maskawa-rotated weak eigenstates (e.g., d' = \Vd,u\^d-\- 
\Vs^u\'^s), normalized such that for v scattering at leading order in QCD, — {l—5ii){d'-\-...). 
The full expressions for jFj appear in, for example, Eq. (10) of Ref. [22]. To include charm 
production, we introduce the parameter A = /{Q^ + Tn1) (i.e. one should simply set A ^ 1 
for the light quark contributions). Our result is that 

/i,3,(e,g^)= '-^^-^f'^'^ (3.24) 

h{i.Q')-^-^f^ (3.25) 

/4(e,Q^)=^^^|^- (3.26) 

To make contact with the measurable structure functions Fj, we instead write our results in 
terms of the partonic functions and convolutions: 

j=l,5 

The quantities Hi and Qi are 

n.{^,Q')= (3.28) 

Gi(^,Q')- f di'Uiii^Q^) (3.29) 

The coefficients a*, /9j and 7j are listed in Tables IV- VI. 

As can be seen from Table IV, the a* coefficients are mostly diagonal, however, there 
is some mixing of the parton functions Ti for F4 and F5. This was already seen in our 
analysis that included target mass effects in the coUinear limit. In fact the terms Y.i oi^j^j 
contributions to Fi in Eq. (3.27) are a factor of 1/(1 + /i^^) suppressed relative to our earlier 
results Eqs. (14)-(18) in Ref. [22] which included target mass corrections in the collinear 
limit. The kx corrections are also responsible for the remaining terms proportional to Tij 
and in Eq. (3.27). 



2. Massless Quark Limit and the NC Process 

Finite charm mass effects are mostly important at intermediate energies. For low energy 
processes, charm production is a negligible effect, whereas at very high energies, charm can 
be treated effectively massless. In these limits, we can evaluate the target mass corrections 
assuming three or four massless quarks. Only three independent helicity amplitudes survive 
the limit of massless quarks, or more generally the limit where the quark fields that build 
the currents in Eq. (2.2) are mass- degenerate. This is discussed in the appendix of Ref. [22]. 
For massless quarks in CC interaction, or generally for the NC interaction we then have: 
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(3 
(3 



30) 
31) 



so /s = /2 and f4 — f2 — fi- Given these relations, one finds that 



W4= -iiWi + -r^W2 



32) 




(3 



33) 



which lead, with the usual definitions, to the expressions 




1 F2 

2 2x 



(3 



34) 



(3 



35) 



where all structure functions are understood to be including the target mass corrections of 
Eq. (3.17). Thus, in the massless quark limit for CC interactions, or for the NC process, the 
evaluation of Fi, F2 and F3 is sufficient to fully describe the target mass corrections, even 
when lepton masses cannot be neglected compared to the incident neutrino energy. 



We have already commented on the comparison between the ^ scaling corrections derived 
from the OPE and those obtained in the collinear parton model in Section HID. In this 
section we parallel this discussion with a numerical evaluation of the charged current neutrino 
scattering structure functions F2 and F3 in both of these approaches. 

Figs. 1-4 show results of F2 and -F3 at = 1 and 4 GeV^. The figures 1-4 (a) show 
the absolute scale of the structure functions, while Figs. 1-4 (6) show ratios of structure 
functions. With (a) and (6) together, one can avoid overstressing kinematic regions with 
large corrections which are, nevertheless, not relevant because the value of the structure 
function is small. The structure functions in the figures are evaluated using the GRV parton 
distribution functions [33], which have 3 active flavors. Here and below, we work in the 
convention of three active flavors and in the MS factorization scheme. 

Fig. 1 (a) shows F2 aX. = 1 GeV^ as a function of x. The solid curve is our full result 
of target mass corrected NLO scahng corrections derived from the OPE) F2(^, 1 GeV^). 
The dot-dashed curve is the NLO corrected ^2^°^, 1 GeV^) without TMC, and the dashed 
curve is the LO F2^'^^ . 

In Fig. 1 (6), we show ratios of F2's at = 1 GeV^, all as a function of x. The solid 
line shows the ratio of the LO F2{i, 1 GeV^) evaluated with the full TMC to the collinear 
approximated TMC (TMC-COL) corrected F2. The collinear approximation is discussed in 
Ref. [22] and is summarized by 



IV. APPLICATIONS 



A. Operator Product Expansion versus Collinear Parton Model 
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i=l,...,5 



(4.1) 



The dot-dashed hne is the ratio of the full TMC corrected F2 to the TMC-COL F2, both at 
NLO. One finds that the collinear approximation does reasonably well in representing the 
full TMC structure function F2. 

The dashed line in Fig 1 (6) is the ratio of the NLO TMC F2 to the NLO F2 in the 
M ^ hmit (Ff ^). The TMC are of order 15% at x = 0.4 and very large above x ~ 0.7. 
The region x > 0.7 is where the NLO ^3°^ is less than 10% of its value ai x — 0.1. 

Figs. 2 (a) and 2 (b) show the same quantities for Q'^ — A GeV^. Figs. 3 (a), 3 (b) and 
4 (a), 4 (b) are the corresponding results for F3. 

The comparative results may be summarized as follows 

• As is well known, target mass corrections become important at low Q and they are 
suppressed at low x. 

• For < a; < 0.9, the coUinear approximation deviates from the exact OPE results 
by at most 10% for F2 and at most 20% for F3. A more detailed examination of the 
adequacy of the collinear approximation for phenomenology will be left to future work. 

As briefly mentioned above, the second point may be reason for some optimism (though 

by no means a proof) that simple kinematic rescalings provide the dominant hadron mass 
effects in other processes as well. A crucial observation here is that a rescaling term is 
enhanced by a derivative factor 



^/.x ^/ X ^ dF 
f (f) - Fi^} = 5 ^ 



{^'x) + ...^0i — ]F'{x) (4.2) 



for a steeply falling function F{^). E.g. for a toy function F{x) oc (1 — x)'^ the rescaling for 
X — > 1 is enhanced as 1/(1 — x) compared to a multiphcative [1 + 0{^)]F{x) correction 
term. 



B. The Paschos-Wolfenstein Relation and Related Observables 

In the following, integrated cross sections 

u,u _ JdE^^p ^'^NC.CC ^{Eud) I20 GeV<yg^,^<180 GeV , . 

will refer to flux-averaged integrals with a cut on hadronic energy as in the experimental 
analysis [12]. We use the (anti-) neutrino fluxes in a step function form as shown in Fig 5. 
We will consider the counting experiment observables 

R''' ^ ^ (4.4) 



'CC 

as well as the Paschos-Wolfenstein [35] relation 
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The question to what extent the approximation in Eq. (4.5) holds in the hght of the high 
precision data of Rcf. [12] has been of some debate recently. For an ideally iso-scalar target 
and under the neglect of charm production components one has that R" = 1/2 — sin^ 0w 
at arbitrary order in QCD as long as there are no (s — s){x) component in the nucleon 
strange sea [14]. For definiteness, we will consider the scattering of neutrinos on an Z — 
26, A — 56 iron target and charm production with rric — 1.3 GeV. We will be using 
parton distributions [33,34] from the literature, though, which do not contain a (s — s){x) 
asymmetry. In a somewhat idealized Mellin moment language, NLO corrections have already 
been estimated in Ref. [14] and we refer the reader to this reference^ for more detailed 
background information. Our evaluation of Eqs. (4.4), (4.5) complements Ref. [14] in that 
we do not make the approximations required for a moment analysis and that we mimic the 
experimental setup to some extent by including the neutrino flux and hadronic energy cut 
in Eq. (4.3). Technically, Eq. (4.3) is evaluated as a 5-dimensional Monte Carlo integral: 



over the kinematic variables and two convolution variables arising from NLO and TM cor- 
rection integrals. Eq. (4.6) integrates numerically over plus type distributions in the pertur- 
bative coefficient functions. In a MC integral, this has to be done with care to achieve the 
required accuracy. 

Our results for R"'"" are summarized in Table VII for two sets (GRV [33] and CTEQ6 
[34]) of parton distribution functions, for a LO^^ or NLO evaluation and for the standard 
model and the anomalous value of the Weinberg angle. The CTEQ6 parton distribution 
analysis also provides a master formula - Eq. (3) in [34] - to estimate the PDF related error 
on a physical observable through a set of error PDFs. Our results can be summarized as 
follows 

• R'' is insensitive to the Weinberg angle and sensitive to NLO corrections 

• R'^ is insensitive to NLO corrections within its sensitivity to the Weinberg angle 

• The impact of PDF uncertainties from pre-determined PDF fits is inconclusive, e.g. the 
error estimate of CTEQ6 for R" does not overlap with the evaluation based on GRV. 
This fact is not too surprising, though, given there are other factors such as scheme 
(number of active flavours) uncertainties and that the error estimate has not been 
tailored to match this particular high precision application. 



■i-^ While completing this article the effect of NNLO moments were presented in Ref. [36] . 

^^The LO evaluation also neglects target mass effects. The charm mass is kept no-zero along the 
LO slow rescaling prescription. 




(4.6) 
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Prom these results, one cannot derive a conclusive estimate of the impact of NLO corrections 
on the analysis [12] and further work will be required. For now, we will restrict ourselves to 
playing the game to treat Eq. (4.5) as a would-be identity and to solve it for the Weinberg 
angle: 



{GRV LO, sin^ ew=0.2227} 
{GRV NLO, sin2 ew=0.2227} 



{CTEQ NLO, sin^ ew=0.2227} 



= 0.2192(3) (4.7) 
= 0.2192(2) (4.8) 
= 0.2196(9) ±0.0005(1) . (4.9) 



The numbers in parentheses are the next digit in the numerical evaluation. For the evaluation 
of R~ , in addition to the values for E!"'^ in Table VII one also needs the values for 



r = ^ 

in Table VIII giving 



(4.10) 



1 -r ^ ' 

The difference between the numerical values in Eqs. (4.7) and Eqs. (4.8) reflects the 
impact of a LO or NLO evaluation of the cross sections entering R~ . It is actually beyond 
the numerical precision of our calculation which is about 10""'. The error quoted with the 
NLO evaluation using CTEQ6M refers to the master formula (3) in Ref. [34] and has to 
be understood as explained in detail in this reference. For the observable R~ we find a 
very robust stability under NLO corrections even under the non-ideal conditions of a non- 
isoscalar target and in a world of massive charm quarks and a massive target. We also 
find a similar stability with regards to PDF variations as long as they do not exploit any 
new physical degree of freedom such as isospin violations or an asymmetry (s — s) {x) in the 
strange sea. From theory alone we cannot investigate if the correlations between and 

in the analysis [12] are correctly represented by the combined ratio R~ . The values in 
Eqs. (4.7) - (4.9), nevertheless, support the expectation that the NuTeV anomaly is not a 
technical NLO effect. 



V. SUMMARY 

In oscillation searches, neutrino cross sections relate event rates with neutrino fluxes 

or the absence of events with bounds on oscillation parameters. A recent measurement of 
the Weinberg angle which deviates from the standard model expectation is based on the 
observation of neutrino interactions on iron. Contrary to the photon exchange case, the 
massive boson propagator is non-divergent at vanishing invariant boson mass. Neutrino 
cross sections are, therefore, dominated by perturbative interactions with quarks and gluons 
provided the neutrino energy is high enough. The dominance of deep inelastic interactions 
typically sets in already above a few GeV neutrino energy. Typical values are not high 
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enoTigh. though, to rely on the lowest order parton model picture for anything but a crude 
estimate. Subleading power suppressed terms and logarithmic corrections are unlikely to be 
separable in the data and have to be considered on the same footing in theory. In this article 
we have revisited, summarized and extended the inclusion of perturbative NLO corrections 
and heavy quark as well as lepton and target mass terms into the twist-2 component of 
weak structure functions and perturbative neutrino cross sections. Further corrections not 
considered here are then, by definition, cither of NNLO perturbative accuracy or of higher 
twist. A factor which we have barely discussed here and which does not seem to lend itself to 
a systematic treatment easily is the impact of our imperfect knowledge of parton distribution 
functions. Obviously, this is even more true for higher twist parton correlations. This article 
may be seen as a technical basis and starting point to address some of these points in more 
detail in the future. We have set up a Monte Carlo integration program for weak structure 
functions and integrated cross sections and, for now, provided numerical results for two 
warm-up applications: First, we have investigated the accuracy of target mass corrections 
in the approximate coUinear parton model which we found to be better than 10% for F2 
and better than 20% for F3. This result is interesting for computational efficiency and as an 
expectation for the accuracy of similar approximations to processes where the OPE is not 
applicable. Second, we have looked into observables that are related to the measurement of 
the weak mixing angle in neutrino scattering. Though we confirmed the expected stability 
of the Paschos-Wolfenstein relation, our first results do not allow any final conclusion on 
the impact of corrective terms to an experimental analysis based on a difi^erent cross section 
model. Future work will be done along this direction which will also take PDF uncertainties 
into account. 

NOTE ADDED ON MOMENT ESTIMATES 

After we submitted this article, in addition to Refs. [14, 36] analytic moment estimates 
of NLO corrections to neutrino cross sections were also presented in Ref . [37] . The authors 

point out a different sign for the NLO correction to W than we find. Under the same 
approximations*** and conventions^^^ as used in [37], we reproduce the results in Eqs. (5.14) 
and (5.15) of that paper. We have then taken the approximate moment approach as a 
starting point and examined the implications of the different approximations it makes. First, 
we find that the sign of the NLO correction to E!" is not flipped by the hadronic energy cut 
[12] in Eq. (4.3). The hadronic energy cut was not included in [37]. Some of the other 
approximations made in [37], however, can flip the sign of the correction to and thus 
explain the seeming discrepancy between our calculation and that of [37]. 



***The approximations required for making moment estimates and obtaining energy-independent 
ratios R^'^ are: (i) neglect of evolution, {ii) Mz,w 00 limit of the boson propagator, (in) neglect 
of the hadronic energy cut. 

^^^I.e., to reproduce Ref. [37] we switch to the 4 flavor DIS factorization scheme {rric = 0) adopted 
there, instead of our canonical 3 flavor MS convention [rric 7^ 0). 
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An approximation that has imphcations for the sign of the NLO correction to is the 
moment approximation itself. The exact fi dependence of the parton distribution functions 
and coefficient functions is parametrically of order a'j.{Q'^) hi{Q'^/Q'^) (for average Q"^ = Q^) 
but its neglect overestimates (underestimates) typical Q values at low x (large x) where 
structure functions are larger (smaller) for increasing Q. This leads to an overestimate 
of neutrino cross-sections of about ~ 10% — 15% around E,y — 100 GeV. In addition, 
for fixed parton distributions (or their second moments), there is a factorization scheme 
dependence in the NLO correction and a comparison of different conventions is of limited 
physical significance. Also, the classification of LO and NLO terms is generally not the same 
ior Uf — 3 or Uf — A active parton fiavours. We find that both the Q :^ Q approximation and 
the dependence on the factorization scheme conventions are capable of flipping the sign of 
the NLO correction to . More relevant of course than the volatile sign of this correction 
is the question - posed in [37] and above in Section IV B - of what impact a complete 
and consistent QCD treatment of neutrino cross sections has on the analysis of the data in 
Ref. [12]. More detailed numerical results related to NuTeV observables will be presented 
in a future publication [38] . 
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APPENDIX A: GENERAL DIS CROSS SECTIONS 

The electroweak weights in Eq. (2.8) for a general {vi — ai) leptonic current are given by 

Eiy [{ml + 4m2mi + m| + 2EiMxy)af + vf{ml — 4m2mi + m| + 2EiMxy)] , . 

^ 4M{Ef - mj) 

[4{y - l)Ef + 2MxyEi + (mi + ms)^] + vf [4.{y - l)Ef + 2MxyEi + (mi - ma)^] . . 
" 8{Ef-ml) 
aiViEi [-ml + ml + 2EiMx{y -2)]y 

= AM{Ef - ml) ^^^^ 

_ af [(mi - m2f + 2EiMxy] (mi + ma)^ + -uf (mi - ma)^ [(mi + ma)^ + 2.g;Ma;y] 
" 8M2(£;f - m?)x ^^^^ 

_ -g; [(mi + m2)(m2 + mi(y - l))af + (mi - m2)(ymi - mi - ma)] . 
~ 2M{Ef - ml) ^^^^ 

For CC neutrino scattering (vi — ai — 1, mi — 0) this reduces to 

""■ = 2111+-^ 

ml Mx\ 

'"^=i-4-(^+2^)^ (^^' 

m^ mnW , . , 
-^=4^ + 2^ 

^5 = _ (AlO) 

^ MEi ^ ' 

In the above Eqs. (A6)-(A10) one has ma = m^ for tau neutrinos and for electron or muon 
neutrinos ma = ni^. ^^ ~ 0. For NC interactions (mi = ma = m) we have 

Eiy [m\?,af - v^) + EiMjaj + )xy] 

= 2(i?f-m^)M 

[2E;(y - 1) + a;|/M] + [2m^ + 2£^f (y - 1) + EiMxy] 

= 4(Ef - m^) ^^''^ 

viaiEfx{y - 2)y 

= - 2(£;f-m^) ^^''^ 
afEim^y / » . 

--- M(W-m^) ^^^'^ 

"^^- M(^f-m^) 

where m = for neutrino neutral current scattering. For charged lepton electromagnetic 
interactions the Wi in Eqs. (All)-(A15) with {a;,Wi>a} — 0, vi — 1 have to be inserted into 
Eq. (2.10). 
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TABLES 







TABLE L Coefficients A} 


in Eq. (3.17). 
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TABLE IL Coefficients B] 


in Eq. (3.17). 
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TABLE in. Coefficients C. 


in Eq. (3.17). 
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TABLE IV. Coefficients a\ in Eq. (3.27). 



a) 


z = 1 


i = 2 


z = 3 


i = 4 


i = i 


X 











J = 2 













J = 3 








2^ 





j = 4 










1 


Xp^(i+ij4') 

2/ix2 




j = 5 












i = 5 




X 



TABLE V. Coefficients /35 in Eq. (3.27). 



/3] i = l i = 2 i = 3 i = A i = 5 



j = l ^ 

j = 2 

i = 3 ^ 

j = 4 ^(2-;.0 '-f 

j = 5 4Mx2a-^x) 0^ 



TABLE VL Coefficients -/j in Eq. (3.27). 
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PDF (sin2 Gw) 






GRV NLO (0.2227) 


0.3120 (0.3120) 


0.3844 (0.3845) 


GRV LO (0.2227) 


0.3125 


0.3860 


GRV NLO (0.2277) 


0.3088 


0.3839 


CTEQ6 NLO (0.2227) 


0.3105 ± 0.0006 


0.3841 ± 0.0038 



TABLE VIL The ratios i?^'^ as defined in Eq. (4.4) eval- 
uated for different parton distributions [33,34], at leading 
and higher order and for two values of the Weinberg an- 
gle. The error quoted for the CTEQ6 PDF refers to the 
master formula (3) in [34] and must be understood as ex- 
plained in detail in this reference. In the first line, the 
numbers in parentheses refer to a pcrturbative expansion 
of the ratios R^'^ directly (instead of the ratios of pertur- 
batively expanded cross sections in Eq. (4.4)); i.e. schemat- 
ically By = R^QS + asi?Q>, inside the parentheses. 



PDF 


r 


GRV NLO 


0.3009 


GRV LO 


0.3014 


CTEQ6 NLO 


0.2909 



TABLE VIII. The ratio r as defined in Eq. (4.10) eval- 
uated for different parton distributions [33,34], at leading 
and higher order. 
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FIG. 1. (a) The charged current neutrino structure function F2 at = 1 GeV^ evaluated 
with NLO ^-scahng corrections (sohd hne) and in LO (dashed) and NLO (dot-dashed) under the 
neglect of target mass corrections. (6) Ratio of the LO ^-scaling evaluation of F2 as in (a) with 

= 1 GeV^ to the LO evaluation in the coUinear approximation is shown by the solid line. The 
dot-dashed line shows the same ratio at NLO. The dashed line shows the ratio of the NLO target 
mass corrected F2 to the NLO F2 in the M — limit. 
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FIG. 3. Same as Fig. 1 for the structure function xF^. 
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FIG. 5. The neutrino (solid) and anti-neutrino (dashed) flux used in Eq. (4.3). 
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